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Abstract 


Background:  Bacterial  infections  of  wounds  impair  healing  and  worsen  scarring.  We  hypothesized  that 
transcriptome  analysis  of  wounds  infected  with  Klebsidla  pneumoniae  (Kp.)  or  Pseudomonas  aeruginosa  {Pa.)  would 
indicate  host  responses  associated  with  the  worse  healing  of  Pa.  than  Kp.  infected  wounds. 

Methods:  Wounds  created  on  post  operative  day  (POD)  0  were  infected  during  the  inflammatory  phase  of  healing 
on  POD3  and  were  harvested  on  POD4  for  microarray  and  transcriptome  analysis.  Other  wounds  received  topical 
antibiotic  after  infection  for  24  hours  to  promote  biofilm  development,  and  were  harvested  on  POD6  or  POD12. 

Results:  Wounds  infected  for  24  hours,  relative  to  uninfected  wounds,  elevated  transcripts  of  immune  response 
functions  characteristic  of  infiltrating  leukocytes.  But  P.a.  infected  wounds  elevated  many  more  transcripts  and 
to  higher  levels  than  Kp.  infected  wounds.  Coincidently,  suppressed  transcripts  of  both  wounds  enriched  into 
stress  response  pathways,  including  EIF2  signaling;  however,  this  was  more  extensive  for  Pa i  infected  wounds, 
including  manyfold  more  transcripts  enriching  in  the  Tell  death'  annotation,  suggesting  resident  cutaneous  cell 
toxicity  in  response  to  a  more  damaging  Pa.  inflammatory  milieu.  The  PODS  wounds  were  colonized  with  biofilm 
but  expressed  magnitudes  fewer  immune  response  transcripts  with  no  stress  response  enrichments.  However, 
elevated  transcripts  of  Pa.  infected  wounds  were  inferred  to  be  regulated  by  type  I  interferons,  similar  to  a  network 
unique  to  P.a.  infected  wounds  on  POD4.  On  POD12,  transcripts  that  were  more  elevated  in  Kp.  infected  wounds 
suggested  healing,  while  transcripts  more  elevated  in  P.a.  infected  wounds  indicated  inflammation. 

Conclusions:  An  extensive  inflammatory  response  of  wounds  was  evident  from  upregulated  transcripts  24  hours 
after  infection  with  either  bacterium,  but  the  response  was  more  intense  for  P.a.  than  K.p.  infected  wounds. 
Coincidently,  more  extensive  down  regulated  transcripts  of  P.a.  infected  wounds  indicated  a  stronger  "integrated 
stress  response"  to  the  inflammatory  milieu  that  tipped  more  toward  cutaneous  cell  death.  Unique  to  P.a.  infected 
wounds  on  POD4  and  PODS  were  networks  inferred  to  be  regulated  by  interferons,  which  may  result  from  intracellular 
replication  of  P.a.  These  data  point  to  specific  downregulated  transcripts  of  cells  resident  to  the  wound  as  well  as 
upregulated  transcripts  characteristic  of  infiltrating  leukocytes  that  could  be  useful  markers  of  poorly  healing  wounds 
and  indicators  of  wound  specific  treatments  for  improving  outcomes. 
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Background 

Wound  healing  proceeds  in  an  orchestrated,  highly  ordered, 
and  overlapping  cascade  of  hemostasis,  inflammation,  prolif¬ 
eration,  and  remodeling  [1,2]  that  can  be  deranged  by  local 
or  systemic  pathologies  such  as  stress  [3,4],  diabetes  [5-9], 
and  infections  [10-13],  including  biofilm  infections  [14,15], 
which  can  result  in  delayed  wound  closure  and  worsened 
scarring. 

Microbial  colonization  of  both  acute  and  chronic 
wounds  is  inevitable.  In  many  settings,  gram-positive  en¬ 
dogenous  skin  flora  predominate  [16].  While  this  is  also 
the  case  for  combat  wounds  [17],  combat  casualties  hos¬ 
pitalized  for  definitive  tertiary  care  in  medical  treatment 
facilities  are  at  high  risk  for  nosocomial  infections 
[18,19]  that  develop  days  after  injury  and  are  largely 
due  to  multi-drug  resistant  gram-negative  organisms  in¬ 
cluding  Acinetobacter,  Pseudomonas ,  Enterobacter,  and 
Klebsiella  [18,20-22]. 

Biofilms  are  highly  differentiated  and  spatially  orga¬ 
nized  three-dimensional  structures  consisting  of  matrix- 
enclosed  communities  of  one  or  more  microbial  species 
that  form  on  colonizable  surfaces  [23].  Biofilms  form 
rapidly  in  acute  wounds  [24,25]  and  contribute  to  the 
pathogenesis  of  chronic  non-healing  wounds  [26-28]. 

Inflammation  is  the  normal  response  to  tissue  injury 
and  microbial  infection,  and  in  proper  degree  is  required 
for  tissue  repair  and  healing.  However,  persistent  inflam¬ 
mation  or  uncontrolled  inflammation  such  as  observed 
in  chronic  and/or  infected  wounds  results  in  poor 
wound  outcomes  due  to  delayed  wound  closure  and 
hypertrophic  scarring  [29].  The  extent  of  inflammation 
and  its  clearance  in  a  timely  fashion  can  influence  the 
rate  and  quality  of  wound  healing  [30,31]. 

Differences  in  the  host-response  to  pathogens  have  been 
associated  with  virulence  mechanisms.  P.  aeruginosa  is  a 
common  gram-negative  bacterium  of  nosocomial  and  life- 
threatening  infections  of  immuno-compromised  patients 
[32].  It  possesses  many  virulence  factors  such  as  exopro¬ 
teases  (e.g.,  elastase),  siderophores,  exotoxins,  hydrogen 
cyanide,  and  pyocyanin  to  attack  host  defenses.  K  pneu¬ 
moniae  account  for  a  substantial  percentage  of  nosocomial 
infections  in  neonates,  patients  undergoing  respiratory 
therapies,  and  patients  hospitalized  in  urology  and  burn 
wards,  but  less  is  known  about  its  virulence  in  wound  in¬ 
fection  and  healing.  Factors  such  as  MagA  (responsible  for 
the  mucoviscosity  that  correlates  with  high  serum  resist¬ 
ance  in  vivo  [33]),  siderophore  aerobactin  [34],  and  cell 
surface-associated  fimbriae  [35]  have  been  implicated  in 
K  pneumoniae  virulence. 

Using  a  full- thickness  dermal  rabbit  ear  wound  model, 
we  previously  demonstrated  that  bacterial  wound  infec¬ 
tions  transitioned  from  active  planktonic-phase  infections 
to  biofilm-dominant  infections  that  delayed  wound  clos¬ 
ure  (epithelialization)  and  granulation  tissue  in-growth 


[11-13].  The  extent  of  wound  healing  impairment  de¬ 
pended  on  the  infecting  bacterial  species.  Pseudomonas 
aeruginosa  (P.a.)  caused  the  most  healing  impairment, 
Klebsiella  pneumoniae  (Kp)  slightly  delayed  closure  rela¬ 
tive  to  uninfected  wounds,  and  Staphylococcus  aureus  was 
intermediate  between  P.a .  and  Kp.  [12].  These  infections 
stimulated  the  expression  of  inflammatory  cytokines  such 
as  interleukin  l-(3  (IL-1|3)  and  tumor  necrosis  factor-a 
(TNF-a)  to  levels  proportional  to  the  wound  healing  im¬ 
pairment  [12]. 

Global  gene  expression  analysis  using  microarrays  has 
been  used  to  define  events  in  tissue  repair  [1],  such 
as  those  that  differ  between  healing  and  non-healing 
venous  leg  ulcers  [36],  skin  wounds  and  oral  mucosa 
wounds  (a  “privileged”  tissue  that  heals  rapidly  without 
scar)  [37],  as  well  as  fetal  wounds  which  heal  without 
scar  and  postnatal  wounds  that  scar  [38]. 

Here,  we  have  used  microarray/transcriptome  analyses 
of  wounds  infected  with  P.a.  or  Kp.  to  identify  similar¬ 
ities  and  differences  in  the  wound  responses  to  these 
differently  virulent  bacteria.  The  findings  contribute  to 
understanding  the  molecular  foundations  of  the  im¬ 
paired  healing  caused  by  bacterial  infections. 

Methods 

Animals 

Under  an  approved  protocol  by  the  Animal  Care  and  Use 
Committee  at  Northwestern  University,  adult  female 
New  Zealand  white  rabbits  (3-6  months,  ~3  kg)  were 
acclimated  to  standard  housing  and  fed  ad  libitum.  All 
animals  were  housed  in  individual  cages  under  constant 
temperature  and  humidity  with  a  12-hour  light-dark  cycle. 

Bacterial  strains  and  culture 

P.a.  strain  PAOl  was  obtained  from  the  laboratory  of 
Dr.  Barbara  H.  Iglewski,  University  of  Rochester  Medical 
Center.  Kp.  strain  BAMC  07-18  was  kindly  provided  by 
LTC  Clinton  Murray  of  Brooke  Army  Medical  Center, 
Fort  Sam  Houston.  To  prepare  the  bacteria  for  the 
wound  infection,  P.a.  PAOl  and  Kp.  BAMC  07-18  were 
grown  in  Luria  (LB)  broth  and  tryptic  soy  broth  (TSB), 
respectively.  After  overnight  incubation  at  37°C,  0.5-mL 
of  the  P.a.  and  Kp.  cultures  were  transferred  to  9.5-mL 
of  fresh  sterile  LB  and  TSB,  respectively,  and  incubated 
at  37°C  until  the  cultures  reached  the  log  phase.  Bacteria 
were  harvested  and  washed  in  phosphate -buffered  saline 
(PBS)  once  by  centrifugation  at  4,000  rpm  for  15-minutes 
at  4°C.  The  resultant  pellets  were  resuspended  in  PBS  and 
adjusted  to  an  optical  density  of  0.5  at  the  wavelength 
of  600-nm  (OD600).  An  OD600  0.5  was  equivalent  to 
108  CFU  per  mL,  which  was  determined  empirically  for 
each  strain  of  bacteria  used.  Ten-pL  of  bacteria  suspen¬ 
sion  at  OD600  0.5  from  each  strain  was  used  as  the 
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inoculum  (approximately  1  x  106  CFU  per  wound)  for  in¬ 
fecting  the  rabbit  ear  wounds. 

Wounding  and  infection 

Wounding,  bacterial  infection,  and  biofilm  formation 
were  adopted  from  our  previously  published  wound  bio¬ 
film  model  [11-13].  Briefly,  rabbits  were  anesthetized 
with  an  intramuscular  injection  of  ketamine  (22.5  mg/ 
kg)  and  xylazine  (3.5  mg/kg)  mixture  prior  to  surgery. 
Ears  were  shaved,  sterilized  with  70%  ethanol,  and 
injected  intradermally  with  1%  lidocaine  /l:  100,000 
epinephrine  at  the  wound  sites.  A  total  of  18  rabbits 
were  used.  Four,  6-mm  diameter,  full-thickness  dermal 
wounds  were  created  in  the  central  area  of  the  ventral 
ear  down  to  perichondrium  (4  wounds  per  ear).  To 
minimize  variations  in  dermal  thickness  and  blood  flow 
to  the  wounds,  distal  and  proximal  areas  of  the  ear 
were  not  used.  After  surgery,  wounds  were  dressed  with 
Tegaderm  (3  M  Health  Care,  St.  Paul,  MN),  a  semi¬ 
occlusive  transparent  film.  Individual  wounds  were  inoc¬ 
ulated  with  106  bacteria  on  postoperative  day  (POD)  3 
and  redressed  with  Tegaderm.  Bacteria  were  allowed  to 
proliferate  in  vivo  under  the  Tegaderm  dressing.  A  top¬ 
ical  antibiotic,  Ciloxan  ointment  (Ciprofloxacin  0.3%, 
Alcon,  Fort  Worth,  TX),  was  applied  on  POD4  to  elim¬ 
inate  free-floating,  planktonic-phase  bacteria,  leaving  a 
predominate  biofilm.  To  prevent  seroma  formation  and 
re-growth  of  planktonic  bacteria,  an  antimicrobial,  ab¬ 
sorbent  dressing  containing  polyhexamethylene  biguan- 
ide  (Telfa  AMD,  Tyco  Healthcare  Group,  Mansfield, 
MA)  was  applied  to  wounds  on  PODs  5,  6,  and  then 
every  other  day  until  harvest.  Uninfected  and  infected 
wounds  were  equally  treated  with  the  antibiotics  and 
dressings.  All  dressings  were  examined  daily  throughout 
the  protocol. 

Harvesting  of  wounds 

After  euthanizing  the  animals  by  intracardiac  euthasol 
injection,  wounds  were  harvested  for  RNA  extraction. 
All  wounds  were  excised  using  a  7-mm  biopsy  punch 
(Acuderm  inc.,  Fort  Lauderdale,  FL).  Wounds  from  dif¬ 
ferent  study  groups  were  harvested  at  POD  4  (prior  to 
the  application  of  local  antibiotic),  POD6,  and  POD  12, 
and  RNA  was  isolated  for  microarray  and  reverse  tran¬ 
scription  quantitative  PCR  (RT-qPCR)  analyses. 

RNA  extraction 

For  mRNA  extraction,  the  dermal  layer  on  the  dorsal 
side  of  the  ear  wound  was  removed  and  the  wound  bed 
was  punched  out  and  samples  of  all  wounds  were  frozen 
in  RNAlater  (Ambion,  Austin,  TX)  at  -80°C  for  molecu¬ 
lar  analysis.  Wound  samples  were  homogenized  using  a 
Mini-bead  beater-8  (Biospec  Products  Inc,  Bartlesville, 
OK)  using  Zirconia  beads  (2.0  mm  diameter,  Biospec 


Products  Inc)  in  the  presence  of  Trizol  Reagent  (Sigma- 
Aldrich,  St.  Louis,  MO).  Total  RNA  was  isolated  according 
to  the  manufacturer  s  protocol.  Contaminating  genomic 
DNA  during  RNA  preparation  was  removed  using  the 
Turbo  DNA- free  kit  (Ambion,  Austin,  TX).  In  general, 
1  pg  of  RNA  was  obtained  from  each  wound.  The  quantity 
and  quality  (Absorbance  260/Absorbance  280  ratio,  28 
s/18  s  RNA  ratio,  RNA  integrity  number)  of  the  extracted 
RNA  was  determined  using  a  NanoDrop  2000C  UV-vis 
Spectrophotometer  (Thermo  Scientific,  Wilmington,  DE) 
and  the  Agilent  2100  Bioanalyzer  (Agilent  Technologies, 
Santa  Clara,  CA),  respectively. 

Microarray 

RNA  (1  pg  each)  was  collected  from  5  control  (unin¬ 
fected),  5  K.p. -infected,  and  5  P.a. -infected  wounds  at 
each  time  point  (POD4,  6,  and  12;  45  samples  total). 
Each  RNA  sample  was  labeled  with  the  One- Color 
Quick- Amp  Labeling  Kit  (Agilent  Technologies)  and  hy¬ 
bridized  to  an  Agilent  custom  rabbit  microarray  (Agilent 
Product  #  G2519F)  containing  4  sets  of  44,000  probes 
per  slide.  The  microarray  contains  43,803  rabbit  probes 
representing  a  total  of  19,938  unique  rabbit  sequences. 
All  labeling  and  hybridization  was  performed  at  the  Uni¬ 
versity  of  Florida  Interdisciplinary  Center  for  Biotechnol¬ 
ogy  Research  according  to  standard  Agilent  recommended 
protocols. 

Data  collection  and  analysis 

Data  preprocessing  was  performed  using  the  Bioconduc¬ 
tor  Limma  package.  The  background  correction  was  per¬ 
formed  using  Normexp  with  offset  value  50.  Quantile 
normalization  was  then  used  to  correct  for  the  non- 
biological  variations  across  arrays.  Probes  were  filtered 
based  on  the  criteria  that  at  least  39  samples  must  have 
a  signal  level  above  background  and  negative  controls. 
The  signal  intensity  values  were  further  normalized  by 
subtracting  (in  log2  scale)  the  time  zero  sample  mean 
(i.e.,  POD3,  immediately  after  wounding).  Differentially 
expressed  genes  were  identified  using  routines  imple¬ 
mented  in  the  Limma  package  to  fit  linear  models  to  the 
normalized  expression  values  [39].  The  variance  used  in 
the  t-score  calculation  was  corrected  by  an  empirical 
Bayesian  method  for  better  estimation  under  small  sam¬ 
ple  size.  All  /^-values  were  corrected  by  the  Benjamini- 
Hochberg  procedure. 

Ingenuity  pathway  analysis 

Ingenuity  Pathway  Analysis  (IPA,  Ingenuity  Systems  Inc.) 
was  used  to  analyze  enrichments  of  DEGs  (infected  vs.  un¬ 
infected  wounds,  and  ICp.-  vs.  P.a. -infected  wounds)  into 
functions,  pathways  and  networks.  Because  IPA  does  not 
analyze  rabbit  gene  identifiers  (IDs),  they  were  first  con¬ 
verted  to  their  human  ortholog  gene  IDs.  Out  of  the 
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19,879  sequences  on  the  microarray,  13,481  unique  se¬ 
quences  mapped  to  human  IDs  and  were  ‘Analysis -Ready 
in  IPA  (68%).  Data  were  analyzed  through  the  use  of 
causal  analysis  approaches  in  IPA  to  infer  upstream  bio¬ 
logical  causes  and  probable  downstream  biological  effects, 
using  enrichment  score  (Fishers  Exact  Test  ^-value)  to 
measure  overlap  between  observed  and  predicted  gene 
sets,  and  z-score  to  match  observed  and  predicted  up/ 
down  regulation  patterns  [40]. 

Reverse-transcription  qPCR 

Subsets  of  RNAs  were  used  for  determining  transcript 
levels  of  specific  target  genes  by  RT-qPCR.  RNA  was 
used  to  prepare  cDNA  using  superscript  II  (Invitrogen) 
with  100-ng  of  random  primers  (Invitrogen).  RT-qPCR 
analyses  using  SYBR  green  was  performed  utilizing  an 
ABI  prism  7000  sequence  detection  system  (Applied 
Biosystems,  Foster  City,  CA).  Genes  that  were  highly 
differentially  regulated  in  either  K.p.  or  P.a.  infected 
wounds  as  shown  in  the  microarray  analyses  were  selected 
for  validation  by  RT-qPCR.  PCR  primers  were  designed 
using  the  Primer  3  program  (http://frodo.wi.mit.edu/).  Ex¬ 
pression  of  each  gene  was  normalized  to  the  level  of  the 
house  keeping  gene,  glyceraldehyde  3-phosphate  dehydro¬ 
genase  (GAPDH),  to  get  ACt.  The  2_AACt  method  was  used 
to  calculate  gene  expression. 

Results 

To  identify  rabbit  ear  wound  transcriptomic  responses 
to  infections  with  K.  pneumoniae  ( K.p .)  or  P.  aeruginosa 
(P.a.),  we  used  the  experimental  design  summarized  in 
Figure  1.  Six-mm  diameter  wounds  were  created  and  in¬ 
fected  3  days  later  (post-operative  day  three,  POD3)  with 
planktonically  grown  K.p.  or  P.a.  (106  colony- forming 
units,  CFU).  At  three  post-infection  time  points,  wound 
RNA  was  harvested  for  microarray  analysis.  Wounds 
were  first  harvested  the  day  after  infection  on  POD4. 
Other  wounds  infected  for  24  hours  were  treated  with 
local  antibiotic  on  POD4  to  kill  off  planktonic  bacteria. 
The  RNA  was  isolated  from  these  biofilm-infected 


wounds  [12]  two  days  later  on  POD6,  and  on  POD  12, 
when  uninfected  wounds  were  fully  epithelialized,  K.p .- 
infected  wounds  were  slightly  open  (~1.7  mm  epithelial 
gap  [12]),  and  P.a. -infected  wounds  were  largely  open 
(-4  mm  epithelial  gap  [12]).  Gene  expression  levels  of 
the  infected  wounds  at  each  time  point  were  compared 
with  uninfected  wounds,  or  the  two  types  of  infected 
wounds  were  compared  each  other,  and  the  differentially 
expressed  transcripts  (DEGs)  were  analyzed  for  associa¬ 
tions  of  with  pathways,  networks,  and  biological  func¬ 
tions  and  disorders  using  the  Ingenuity  Pathway  Analysis 
(IPA)  database  and  software. 

Tallies  of  differentially  expressed  genes  (DEGs)  in  infected 
vs.  uninfected  wounds 

We  found  1,378  differentially  expressed  sequences  (fold 
change  >  2  with  adjusted  j?-value  <  0.05)  in  infected  as 
compared  to  uninfected  (I/U)  wounds  out  of  the  19,879 
unique  sequences  on  the  microarray;  i.e.,  7%  of  all  gene 
probes  were  I/U  differentially  expressed  (I/U-DEGs)  by 
>2-fold,  when  tallying  across  all  PODs  for  both  P.a.-  and 
K.p.-  infected  wounds  (Table  1). 

Averaging  over  all  PODs,  P.a.-  infected  wounds  con¬ 
tained  3-fold  more  I/U-DEGs  than  K.p.- infected  wounds. 
This  was  largely  due  to  the  5-fold  more  I/U-DEGs  in 
P.a.-  than  TQz -infected  wounds  on  POD4  the  day  after 
bacteria  were  inoculated  into  the  wounds  and  before 
antibiotic  treatment  [12],  At  this  time  exudate  was  abun¬ 
dant  and  the  number  of  >2-fold-changed  I/U-DEGs  was 
219  for  K.p.  -infected  wounds  and  1,113  for  P.a. -infected 
wounds  (Table  1).  There  were  many  fewer  DEGs  on 
POD6  (K.p.:  17;  P.a.:  30)  following  two  days  of  antibiotic 
treatment  when  -1  to  2  logs  fewer  CFU  were  recovered 
than  were  inoculated  (our  previous  study  [12]).  On 
POD  12,  the  I/U-DEGs  had  increased  (K.p.:  162;  P.a.: 
161,  Table  1). 

The  complete  list  of  all  19,879  gene  probes  and  associ¬ 
ated  normalized  log2  fold-changes,  adjusted  ^-values, 
and  updated  annotations  is  provided  in  Additional  file  1, 
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Table  1  Summary  of  genes  with  >2-fold  differential 
expression  in  K.p.-  or  P.a-infected  wounds  vs.  uninfected 
wounds 


POD4 

POD6 

POD12 

Infecting  bacteria: 

K.p. 

P.a. 

K.p. 

P.a. 

K.p. 

P.a. 

Upregulated  genes*: 

168a 

701 

12 

21 

77 

129 

Downregulated  genes*: 

51 

412 

5 

9 

85 

32 

Totals 

219 

1113 

17 

30 

162 

161 

"Based  on  >  2  fold  change  of  probe  intensity  with  an  adjusted  p  value  <  0.05. 
""Overlapping  l/U  DEGs  are  included. 


and  these  data  have  been  deposited  in  the  NCBI  GEO 
database  (Accession  number  GSE51167). 

Few  l/U-DEGs  overlapped  between  PODs,  but  many 
l/U-DEGs  were  common  to  K.p.-  and  P.a-infected 
wounds  on  each  POD 

To  compare  overlapping  genes  between  PODs  and  infec¬ 
tion  types  on  each  POD,  we  used  the  top  200  differentially 
expressed  genes,  consisting  of  the  100  most  up-regulated 
(Upr)  and  100  down-regulated  (Dnr)  I/U-DEGs,  in  order 
to  compare  equal  numbers  of  I/U-DEGs.  Not  a  single  gene 
overlapped  between  all  three  PODs,  for  either  species,  and 
less  than  7%  overlapped  between  any  two  PODs,  indicat¬ 
ing  the  largely  different  states  of  infection  and  healing 
of  the  wounds  on  the  three  PODs.  In  contrast,  many 
I/U-DEGs  were  common  to  K.p.-  and  P. a. -infected 
wounds  (Figure  2)  indicating  similar  host-responses  to  the 
different  infections.  To  go  beyond  comparing  similarities 
and  differences  of  single  genes,  we  analyzed  similarities 
and  differences  of  pathways  and  functions  associated  with 
the  I/U-DEGs  common  and  unique  to  the  two  infected 
wound  types. 


POD4:  Upr-l/U-DEGs  common  to  K.p.-  and  P.a-infected 
wounds  enriched  in  inflammatory-response  functions 

The  common  I/U-DEGs  to  K.p.-  and  P.a. -infected 
wounds  enriched  in  the  IPA  Function  Annotation  cat¬ 
egory  Inflammatory  Response'.  Downstream  Effects 
Analysis  of  these  common  genes  identified  Functions 
Annotations  whose  activation  state  was  inferred  to  in¬ 
crease,  based  on  literature  compiled  in  the  Ingenuity® 
Knowledge  Base,  to  include:  adhesion  of  immune  cells', 
synthesis  of  reactive  oxygen  species',  recruitment  of  neu¬ 
trophils',  activation  of  myeloid  cells',  activation  of  leuko¬ 
cytes',  production  of  reactive  oxygen  species',  adhesion  of 
mononuclear  leukocytes!  activation  of  phagocytes',  ‘mi¬ 
gration  of  neutrophils',  ‘cell  movement  of  myeloid  cells', 
and  ‘inflammatory  response'.  The  Z-scores  for  enrich¬ 
ment  of  the  I/U-DEGs  in  these  functions  were  between 
3.2  and  2.4  (significance  cut-off  2.0),  predicting  an  acti¬ 
vated  state.  The  top  four  inferred  upstream  regulators  for 
these  Upr  I/U-DEGs  common  to  both  infected  wound 
types  were  lipopolysaccharide  (LPS),  TNF,  IFNG  and  IL1B 
(z-scores  3.5  to  2.5).  These  results  show  that  the  Upr  I/U- 
DEGs  common  to  the  two  infected  wound  types  were 
mostly  involved  in  the  inflammatory  response. 

POD4:  99%  of  I/U-DEGs  were  concordantly  Upr  or  Dnr  in 
K.p.-  and  P.a-infected  wounds 

Seeking  to  identify  differences  between  the  host-responses 
to  the  two  infected  wound  types,  we  mined  more  deeply 
into  the  dataset,  analyzing  all  I/U-DEGs  on  POD4  with  a 
/7-value  of  <0.05  (Benjamini-Hochberg).  The  I/U-DEGs 
common  to  both  K.p.-  and  P.a. -infected  wounds  summed 
to  907  (i.e.,  all  overlapping  I/U-DEGs).  All  of  these  I/U- 
DEGs,  except  for  1%,  or  9,  showed  the  same  direction  of 
differential  expression  in  both  infected  wounds  relative  to 
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uninfected  wounds  (45%  Upr  and  54%  Dnr).  This  result 
suggests  biological  significance  of  even  small  I/U  fold- 
changes  and  indicates  largely  similar  responses  of  the 
wounds  to  the  two  types  of  bacterial  infection. 

The  great  majority  of  Upr-transcripts,  91%,  were  up- 
regulated  to  a  greater  extent  in  P.a.- than  K.p. -infected 
wounds.  Many  of  these  are  mainly  expressed  in  leuko¬ 
cytes,  for  example:  the  super  oxide-producing  NAD  PH 
oxidase  expressed  in  neutrophils,  NCF1,  2  and  4  (2.8-, 
3.8-  and  2.9-fold);  the  G-protein  coupled  receptor  EMR2 
(2.9-fold);  S100A9  (1.9-fold);  CSF3R  (colony  stimulating 
factor  3  receptor,  granulocyte  (2.9-fold);  CSF2,  colony 
stimulating  factor  2,  granulocyte-macrophage  (5-fold);  and 
CSF3,  colony  stimulating  factor  3,  granulocyte  (3.6-fold). 
Conceptually,  as  compared  to  K.p. -infected  wounds,  some 
elevated  transcripts  of  P.a. -infected  wounds  could  come 
from  greater  numbers  of  infiltrating  leukocytes.  Similar 
numbers  of  leukocytes  or  cutaneous  cells  could  have  up- 
regulated  transcription  to  a  greater  extent  in  the  absence  of 
greater  immune  cell  infiltration,  but  this  seems  less  likely. 

The  majority  of  Dnr  I/U-DEGs,  73%,  were  also  down- 
regulated  to  a  greater  extent  in  P.a. -infected  wounds. 
More  Dnr  genes  and  a  greater  extent  of  downregulation 
of  most  overlapping  I/U-DEGs  in  P.a. -infected  wounds 
could  represent  a  stronger  stress  response  of  cutaneous 
cells  to  greater  inflammation  caused  by  P.a.  infection. 

POD4:  Transcripts  upregulated  to  a  greater  extent  in 
K.p.- infected  wounds  appear  to  come  from  cutaneous 
cells  responding  to  the  inflammatory  milieu 

Of  the  9%,  or  36,  of  the  concordantly  Upr-transcripts 
that  were  upregulated  to  a  greater  extent  in  K.p. -infected 
wounds,  only  KRT8,  CYR61  and  CD5  were  enriched  in 
the  function  ‘inflammatory  response.  However,  KRT8 
and  CYR61  are  likely  not  from  leukocytes.  KRT8  (kera¬ 
tin  8)  is  a  type  II  intermediate  filament  that  forms 
intermediate-sized  filaments  in  the  cytoplasm  of  epithe¬ 
lial  cells  [39].  And  in  fibroblasts,  CYR61  appears  to  up- 
regulate  the  expression  of  a  number  of  genes  involved  in 
angiogenesis,  inflammation,  and  matrix  remodeling  dur¬ 
ing  wound  healing  [39].  The  top  ‘Physiological  System 
Development  and  Function  was  ‘Connective  Tissue  Devel¬ 
opment  and  Function  with  seven  enriching  genes  (CYR61, 
IQGAP1,  NF2,  KRT8,  ALDH1A3,  MXD1,  and  CYP27B1). 
Also  another  keratin,  KRT78,  was  among  the  36  genes. 
These  enrichments  suggest  the  genes  more  upregulated  in 
K.p.- infected  wounds  came  from  cutaneous  cells. 

The  second  topmost  ‘Molecular  and  Cellular  Function' 
into  which  the  36  Upr  genes  enriched  was  ‘Cell  Death 
and  Survival'.  Eleven  genes  enriched  in  this  function 
(CD5,  AIFM3,  CYR61,  EGLN3,  ALDH1A3,  KRT8, 
MTFP1,  MXD1,  NF2,  PDK1,  and  SLC5A8).  The  topmost 
Upr  genes  to  a  greater  extent  in  K.p.-  than  P. a. -infected 
wounds  were  CYR61,  EGLN3,  and  HILPDA.  These  three 


genes  are  all  hypoxia  inducible.  Additionally,  the  top  ca¬ 
nonical  pathway  enrichment  of  the  36  genes  was  glycolysis. 

The  36  genes  that  are  Upr  to  a  greater  extent  in  Kp.- 
infected  wounds  are  listed  in  Additional  file  2.  Their  lack 
of  extensive  enrichment  in  inflammation  functions  and 
their  enrichment  in  functions  related  to  surviving  the 
inflammatory  environment  suggest  they  represent  the 
cutaneous  cell  response  to  the  inflammatory  milieu.  Up¬ 
stream  regulators  with  a  predicted  activated  state  for 
regulating  these  36  genes  were  only  two:  HIF1A  and 
NUPR1.  HIF1A  is  the  master  regulator  of  the  response 
to  hypoxia;  under  hypoxia  it  activates  transcription  of 
over  40  genes  [41].  NUPR1  is  as  a  component  of  the 
NUPR1/RELB/IER3  survival  pathway  that  binds  to  chro¬ 
matin  to  change  gene  expression  that  is  associated  with 
resistance  to  stress  [41]  and  inhibition  of  apoptosis  [42]. 

As  an  additional  analysis,  we  directly  compared  differ¬ 
ential  expression  between  K.p.-  and  P.a. -infected  wounds 
(i.e.,  K/P-DEGs),  as  opposed  to  the  previous  analyses 
of  differential  expression  relative  to  uninfected  wounds 
(i.e.,  I/U-DEGs).  The  K/P-DEGs  that  were  more  elevated 
in  K.p. -infected  wounds  summed  to  836  (p-w alue  <0.05, 
Benjamini-Hochberg).  The  topmost  four  canonical  path¬ 
way  enrichments  of  these  were  ‘Cell  Cycle  Control  of 
Chromosomal  Replication',  ‘Cell  Cycle:  G2/M  DNA 
Damage  Checkpoint  Regulation!  ‘Mitochondrial  Dysfunc¬ 
tion',  and  ‘Role  of  BRCA2  in  DNA  Damage  Response'. 
The  topmost  activated  function  was  cell  survival'  (p-value 
5.08E-03;  89  molecules;  activation  z-score  5.66).  The  top¬ 
most  inhibited  function  was  ‘organismal  death'  (p-w alue 
1.96E-03;  150  molecules;  activation  z-score  -11.38).  Similar 
results  were  obtained  using  the  more  limited  set  of  K/P- 
DEGs  with  the  cutoff  of  >2-fold  change.  None  of  the  func¬ 
tion  enrichments  were  related  to  inflammation  or  response 
to  infection.  One  possible  explanation  for  this  result  is 
greater  numbers  of  leukocytes  had  infiltrated  P. a. -infected 
wounds,  importing  more  inflammation  transcripts  than  in 
K.p. -infected  wounds.  These  data  suggest  that  after  24  hours 
of  active  infection,  cutaneous  cells  of  K.p.  -infected  wounds 
activated  pathways  for  surviving  the  particular  conditions  of 
the  infection  such  as  hypoxia  and  reactive  oxygen  species. 

POD4:  Cutaneous  cells  are  suggested  to  activate  an 
"Integrated  stress  response" 

Since  the  response  to  infection  involves  cells  infiltrating 
the  wound  that  only  contribute  transcripts  to  the  wound, 
Dnr- transcripts  come  from  cutaneous  cells.  The  27%,  or 
133,  I/U-DEGs  that  were  downregulated  to  a  greater  ex¬ 
tent  in  K.p.-  than  P.a. -infected  wounds  had  as  their  top¬ 
most  enrichment  in  the  category  ‘Physiological  System 
Development  and  Function'  the  function  annotation  ‘Hair 
and  Skin  Development  and  Function'  (p-value  6.20E-06 
to  1.95E-02;  7  molecules),  additionally  suggesting  these 
transcripts  come  from  cutaneous  cells.  These  133  I/U- 
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DEGs’  topmost  significant  canonical  pathway  enrichments 
were  ‘EIF2  Signaling!  ‘Regulation  of  eIF4  p70S6I<  Signaling’ 
and  ‘mTOR  Signaling’.  Additionally,  the  topmost  activated 
upstream  regulator  inferred  for  the  133  genes  was  the  group 
N-cor  (Activation  z-score  =  2.0),  consisting  of  NCOR1  and 
NCOR2,  nuclear  receptor  co-repressors  that  mediate  basal 
transcriptional  silencing  [43].  In  sum,  these  results  suggest 
that  cutaneous  cells  of  K.p. -infected  wounds  activated  an 
“integrated  stress  response”  which  involves  downregulation 
of  general  translation  and  upregulation  of  specific  transcripts 
to  survive  the  stressors  of  Kp.  infection  [44,45]. 

The  integrated  stress  response  was  also  suggested  for 
P.a. -infected  wounds  as  the  top  three  canonical  pathway 
enrichments  for  212  unique  Dnr  DEGs  of  P.a. -infec¬ 
ted  wounds  were  ‘Mitochondrial  Dysfunction!  ‘Oxidative 
Phosphorylation’  and  EIF2  signaling.  These  as  well  as 
‘Regulation  of  eIF4  and  p70S6I<  Signaling’  and  ‘mTOR  Sig¬ 
naling’  were  also  the  top  5  canonical  pathways  into  which 
enriched  85  Dnr-I/U-DEGs  common  to  K.p .-  and  P.a- in¬ 
fected  wounds.  Although  these  data  indicate  that  the  cuta¬ 
neous  cells  of  both  types  of  infected  wounds  had  activated 
an  “integrated  stress  response”  to  survive  the  inflamma¬ 
tory  milieu,  the  P.a. -unique  I/U-DEGs  enriched  more 
into  the  ‘Cell  Death’  function  category.  Using  212  unique 
Dnr-I/U-DEGs  to  each  wound  type,  17  from  K.p. -infected 
wounds  and  72  from  P.a.  -infected  wounds  enriched  in 
the  ‘Cell  Death’  function  annotation  category. 

POD4:  Discordantly  regulated  transcripts  between 
K.p.-  and  P.a-infected  wounds. 

Of  the  9  (1%)  discordantly  regulated  I/U-DEGs  between 
the  two  types  of  infected  wounds,  four  were  significantly 
Upr  in  P.a.  and  Dnr  in  K.p.:  SPTAN1,  ATP6V1F, 
ATG4A,  and  EEA1  (K/P-DEGs  were  3.8-,  3-,  1.8-  and 
1.8-fold  higher  in  P.a.,  respectively).  Each  of  these  genes 
has  a  function  related  to  intracellular  organelles  [41]. 
SPTAN1  encodes  the  filamentous  cytoskeletal  protein 
spectrin  (alpha  subunit,  aka  Fodrin),  which  functions 
as  a  membrane-stabilizing  scaffold  required  for  the 
function  of  intracellular  organelles.  ATP6V1F  encodes  a 
component  of  vacuolar  ATPase  (V-ATPase),  which  me¬ 
diates  acidification  of  eukaryotic  intracellular  organelles. 
ATG4A  encodes  a  cysteine  protease  required  for  autoph- 
agy  and  is  implicated  in  the  processing  of  mature  autopha¬ 
gosomes  to  facilitate  their  delivery  to  the  lysosomal 
network.  The  EEA1 -encoded  protein  binds  phospholipid 
vesicles  containing  phosphatidylinositol  3-phosphate  and 
participates  in  endosomal  trafficking.  Notably,  P.a.  viru¬ 
lence  factor  ExoS  (expressed  by  PA01  in  our  in  vitro  ex¬ 
periments,  unpublished  observation)  and  in  particular  its 
ADP-ribosyltransferase  activity,  is  required  for  intravacuo- 
lar  replication  of  P.a.  in  epithelial  cells  [46]. 

The  5  I/U-DEGs  that  were  elevated  in  K.p.-  and  sup¬ 
pressed  in  P.a. -infected  wounds  were  LMOl,  BRMS1L, 


CHD1L  (K/P-DEGs  were  2.5-,  1.4,  and  1.5-fold  higher  in 
K.p.),  as  well  as  two  transcribed  loci  EB374055  and 
DN885016.  The  three  known  genes  encode  chromatin 
binding  proteins  [41].  The  LMOl-encoded  protein  is 
suggested  to  regulate  transcription  by  competitively 
binding  to  specific  DNA-binding  transcription  factors. 
BRMS1L  encodes  a  component  of  histone  deacetylase 
complexes.  CHD1L  encodes  a  DNA  helicase  that  plays  a 
role  in  chromatin-remodeling  following  DNA  damage. 
One  possibility  is  that  K.p.  -infected  wounds  retained 
capacity  to  remodel  chromatin  in  response  to  the  stress 
of  the  infectious  inflammation. 

POD4:  Network  with  upstream  regulation  by  type  I 
interferons  unique  to  P.a-infected  wounds 

Seeking  to  identify  differences  in  upstream  regulators  as¬ 
sociated  with  the  gene  expression  patterns  of  the  infected 
wound  types,  we  analyzed  upstream  activating  regulators 
of  bacteria.  They  were  markedly  similar  between  the  two 
wound  types.  Given  the  extensive  overlap  of  identical  and 
similar  genes,  to  further  search  for  discriminating  features 
of  the  wound-response  to  P.a.  and  K.p.  infections,  we  ana¬ 
lyzed  I/U-DEGs  grouped  by  predicted  activated  upstream 
regulators  (z-score  >  2)  in  the  categories,  ‘Cytokine!  ‘G -pro¬ 
tein-coupled  receptor!  ‘Growth  factor!  ‘Ligand-dependent 
nuclear  receptor!  ‘Transmembrane  receptor!  and  ‘Other’. 
The  I/U-DEGs  of  each  wound  type  without  direct  regula¬ 
tory  connections  to  other  molecules  in  the  same  group, 
based  on  the  IPA  database,  were  excluded  from  the  ana¬ 
lysis.  There  were  89  upstream  regulators  of  the  134  mole¬ 
cules  common  to  both  infected  wound  types,  69  upstream 
regulators  of  the  113  molecules  unique  to  P.a-infected 
wounds,  and  3  predicted  upstream  regulators  of  the  66 
molecules  unique  to  K.p. -infected  wounds. 

Between  the  K.p. -P.a.-  common  and  P.a -unique  I/U-DEGs, 
44  upstream  regulators  overlapped.  The  remaining  upstream 
regulators  summed  to  25  and  are  inferred  to  regulate  56 
downstream  I/U-DEGs  unique  to  P.a-infected  wounds.  Fea¬ 
tures  of  this  network  include  regulation  by  type  I  and  type  II 
interferons,  and  16  downstream  molecular  targets  with  more 
than  four  regulatory  connections:  INFG,  STAT1,  STAT4, 
PIM1,  PIM2,  CCL5,  CSF2,  MCL1,  MMP1,  MMP3,  IRF1, 
IRF7,  RSAD2,  EGR1,  DHX58,  and  OASL  (Additional  file  3). 
Although  predicted  to  be  activated  by  type  I  interferons,  nei¬ 
ther  the  type  I  interferons  (IFNA,  IFNA2,  IFNA5,  IFNL3, 
and  IFNA14)  nor  the  type  I  interferon  receptor,  IFNAR1, 
that  were  on  our  microarray  were  upregulated. 

POD6  wounds  have  fewer  inflammation  and  stress-response 
transcripts,  but  Upr-transcripts  from  P.a-infected  wounds 
form  a  network  with  predicted  upstream  regulation  by  type 
I  interferons  as  on  POD4 

By  POD6,  wounds  infected  for  24  hours  and  then 
treated  with  antibiotic  for  two  days  were  likely  colonized 
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with  biofilm,  as  approximately  105  and  104  CFU  of  K.p. 
and  P.a.,  respectively,  were  recovered  from  the  wounds 
(our  previous  study  [12]),  which  is  1  to  2  logs  fewer  than 
the  106  CFU  inoculated  on  POD3.  Although  about  10- 
fold  fewer  P.a.  CFUs  were  recovered  per  wound,  P.a.  is 
known  to  replicate  in  vacuoles  of  epithelial  cells  [46]. 
Consistent  with  low  bacterial  burden  as  biofilm,  both  in¬ 
fected  wound  types  had  far  fewer  I/U-DEGs  on  POD6 
than  on  POD4  (Table  1). 

To  assess  similarities  and  differences  in  the  response  of 
the  wounds  to  colonization  with  K.p.  vs.  P.a.  biofilm,  we 
analyzed  enrichments  of  the  common  and  the  unique  I/U- 
DEGs  of  K.p.-  and  P.a.- infected  wounds.  We  analyzed  the 
top  fold-change  I/U-DEGs,  100  Upr  and  100  Dnr  that 
were  significant  (adjusted  j^-value  <  0.05),  in  order  to  com¬ 
pare  nearly  equal  numbers  of  I/U-DEGs,  since  I/U-DEGs 
with  >2  fold-change  were  so  few.  Neither  K.p.-P.a.-c om- 
mon-I/U-DEGs  nor  the  K.p. -unique  I/U-DEGs  enriched 
in  functions  with  z-score  >2.0.  In  contrast,  the  94 
I/U-DEGs  unique  to  P.a. -infected  wounds  enriched  in 
functions  with  predicted  activation  (z-score  2.7  to  2.0):  ac¬ 
tivation  of  T  lymphocytes!  antiviral  response!  activation  of 
cells!  activation  of  antigen  presenting  cells!  cytotoxicity  of 
lymphocytes!  ‘killing  of  cells!  cell  death  of  immune  cells! 
activation  of  leukocytes!  activation  of  cytotoxic  T  cells! 
antiviral  response  of  cells!  ‘response  of  mononuclear 
leukocytes!  ‘fragmentation  of  DNA!  ‘cytotoxicity  of  T  lym¬ 
phocytes!  quantity  of  cytotoxic  T  cells!  ‘maturation  of 
phagocytes!  ‘maturation  of  antigen  presenting  cells!  ‘im¬ 
mune  response  of  cells!  and  ‘activation  of  macrophages'. 
While  the  147  I/U-DEGs  of  K.p. -infected  wounds  barely 
made  any  regulatory  connection  amongst  themselves,  the 
94  I/U-DEGs  of  P. a. -infected  wounds  connected  into  a 
25-molecule  regulatory  network  in  which  the  major  hubs 
were  STAT1  and  IRF7,  with  each  having  16  direct  regula¬ 
tory  connections.  This  25-molecule  network  had  several 
genes  in  common  with  the  upstream-interferon-predicted 
network  of  I/U-DEGs  unique  to  P.a.- infected  wounds  on 
POD4.  These  common  genes  between  the  POD4  and 
POD6  networks  are  STAT1,  IRF7,  and  DHX58.  Addition¬ 
ally,  the  POD6  network  includes  PSMB10,  GBP1,  and 
OAS2,  while  the  POD4  network  includes  their  highly  re¬ 
lated  molecules  PSMB9,  GBP4,  and  OASL.  These  25  mol¬ 
ecules'  topmost  five  enrichments  into  canonical  pathways 
are  shown  in  Table  2.  The  molecule  names,  and  the  log  ra¬ 
tios  and  p-y alue  of  their  I/U-DEGs  are  listed  in  Additional 
file  4.  Furthermore,  the  Upr  I/U-DEGs  enriched  as  their 
topmost  network:  ‘ Infectious  Disease ,  Inflammatory  Re¬ 
sponse,  and  Antimicrobial  Response  (p-y alue  lxlCT47  right 
tailed  Fisher's  exact  test).  Table  3  lists  the  topmost  five  up¬ 
stream  mammalian  receptors  inferred  to  activate  these 
I/U-DEGs.  Although  type  I  interferons  are  prominent 
among  these,  neither  they  nor  their  receptor  were  upregu- 
lated  in  our  experiments. 


Analysis  of  the  Dnr  I/U-DEGs  for  associations  with 
upstream  regulators  and  downstream  functions  found 
that  despite  the  presence  of  upregulated  inflammatory/ 
antimicrobial  transcripts,  activation  of  stress-response 
pathways  was  not  evident  as  on  POD4. 

POD12:  More  inflammation  and  inhibited  healing 
suggested  for  P.a-infected  wounds 

On  POD  12,  uninfected  wounds  had  completely  epithelia- 
lized,  K.p. -infected  wounds  were  nearly  completely  epithe- 
lialized  (epithelial  gap  ~1.7  mm),  and  P.a.-  infected 
wounds  were  more  open  (epithelial  gap  -4  mm)  [12].  In 
both  K.p.-  and  P.a. -infected  wounds,  biofilm  was  seen  on 
POD12  using  scanning  electron  microscopy  [12]. 

The  I/U-DEGs  common  to  K.p.-  and  P.a.- infected 
wounds  (Upr  and  Dnr)  enriched  in  the  functions  ‘Vasculo- 
genesis'  and  ‘Angiogenesis!  predicted  to  be  activated  with 
z-scores  2.2  and  2.0  (p- value  1.22E-02  and  1.97E-03,  re¬ 
spectively).  The  top-scoring  function  was  ‘Connective  Tis¬ 
sue  Disorders'  (^-values  1.27E-06  to  2.33E-03),  which 
included  collagen  type  XI,  alpha  1!  collagen,  type  XII, 
alpha  1!  collagen,  type  V,  alpha  2! ‘fibrillin  1', ‘interleukin  18 
binding  protein!  ‘lysyl  oxidase-like  2!  ‘matrix  metallopepti- 
dase  1  (interstitial  collagenase)!  ‘thrombospondin  2'.  All 
were  downregulated  except  IL18BP  and  MMP1,  which 
were  upregulated. 

We  also  analyzed  K.p.-  versus  P.a. -infected  wound  dif¬ 
ferential  expression  (K/P).  Transcripts  significantly  ele¬ 
vated  in  P.a.-  relative  to  TQz -infected  wounds  summed 
to  118.  These  molecules  significantly  enriched  into  ‘In¬ 
flammatory  Response',  ‘Connective  Tissue  Disorders',  ‘In¬ 
flammatory  Disease',  ‘Cellular  Movement',  ‘Hematological 
System  Development  and  Function'  and  ‘Immune  Cell 
Trafficking'  (movement  and  migration  of  phagocytes) 
as  topmost  functions.  These  118  genes  are  shown  in 
Additional  file  5.  Notably,  the  gene  with  the  largest  fold- 
change  was  PLI<5,  which  was  10-fold  elevated  in  P.a.-  as 
compared  to  K.p. -infected  wounds.  PLI<5  was  upregu¬ 
lated  in  serum-starved  fibroblasts  as  they  exited  the 
cell  cycle,  and  PLI<5  overexpression  arrested  cells  at  a 
GO/Gl-like  stage  [47].  The  topmost  predicted  upstream 
regulator  for  the  118  molecules  elevated  in  P. a. -infected 
wounds  was  FGF1  (fibroblast  growth  factor  1,  acidic), 
which  trended  toward  an  inhibited  activation  state 
(z-score  -0.818;  /?-value  3.85E-04). 

Conversely,  molecules  elevated  in  K.p.-  relative  to  P.a .- 
infected  wounds  summed  to  86  (Additional  file  6).  A 
single  activating  upstream  regulator  was  the  synthetic 
glucocorticoid  methylprednisone  (z-score  2.45,  j^-value 
8.88E-03,  7  genes:  SCAMPI,  LDLR,  HSD17B12,  SDC2, 
CSNK2B,  ATIC  and  FAH).  Notably,  endogenous  gluco¬ 
corticoid  has  been  suggested  to  suppress  granulation  tis¬ 
sue  formation,  as  mice  with  a  DNA-binding-defective 
glucocorticoid  receptor  had  early  wounds  with  enlarged 
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Table  2  Top  five  canonical  pathway  enrichments  of  25  molecules  from  the  network  on  POD6  unique  to  P.a.-infected 
wounds 


Name 

p  value 

Ratio 

Activation  of  IRF  by  Cytosolic  Pattern  Recognition  Receptors 

1.82E  08 

5/72  (0.069) 

Crosstalk  between  Dendritic  Cells  and  Natural  Killer  Cells 

1.1 1 E  07 

5/96  (0.052) 

Role  of  Pattern  Recognition  Receptors  in  Recognition  of  Bacteria  and  Viruses 

1.54E  07 

5/106  (0.047) 

Role  of  RIG1  like  Receptors  in  Antiviral  Innate  Immunity 

3.22E  07 

4/49  (0.082) 

Communication  between  Innate  and  Adaptive  Immune  Cells 

6.08E  06 

4/109  (0037) 

granulation  tissue  and  a  high  fibroblast  density  [48]. 
PDGF  BB  was  also  a  predicted  upstream  regulator, 
though  it  did  not  reach  the  activating  z-score  cutoff  of 
2.0  (z-score  1.41,  p-vzhie  1.22E-02,  5  genes:  SLC7A7, 
MGST2,  THBD,  FMOl  and  LDLR).  Additionally,  EVPL 
(Envoplakin)  was  2.8-fold  upregulated  in  K.p. -infected 
wounds  over  uninfected  wounds  and  was  2-fold  more 
highly  expressed  in  ICp.-  than  P.a. -infected  wounds. 
Envoplakin  is  a  component  of  the  keratinocyte  cornified 
envelope,  and  is  suggested  to  link  it  to  desmosomes  and 
intermediate  filaments  [39]. 

q-RT-PCR  validation  of  microarray  results 

Analysis  of  gene  expression  using  q-RT-PCR  of  the  same 
RNA  samples  as  used  for  microarray  analysis  validated 
the  microarray  results  (Additional  file  7:  Table  SI).  The 
selected  genes  represented  those  that  were  highly  regu¬ 
lated  either  in  K.p .  or  P.a .  infected  wounds.  The  qRT- 
PCR  results  confirmed  the  trend  in  expression  of  these 
selected  genes  revealed  by  microarray. 

Discussion 

Wound  gene  expression  patterns  differed  greatly  be¬ 
tween  PODs  4,  6  and  12.  This  was  expected  as  the 
wound  states  are  different  on  these  days,  and  respect¬ 
ively  involve  acute  bacterial-induced  inflammation, 
stabilization  of  biofilm,  and  healing.  However,  the  wound 
responses  to  K.p.  and  P.a.  infections  showed  similarities, 
especially  on  POD4.  But  there  were  significant  differ¬ 
ences.  On  POD4,  the  Upr  gene  expression  patterns 


indicating  inflammation  were  largely  similar  between  the 
two  infected  wound  types.  But  P.a. -infection  induced  a 
more  intense  response,  with  Dnr-transcripts  indicating  a 
stronger  stress-response  and  the  death  of  cutaneous 
cells.  The  greater  inflammatory  response  to  P.a.  infec¬ 
tion  is  likely  associated  with  its  pathogenic  virulence  that 
could  have  resulted  in  a  larger  bacterial  load  as  com¬ 
pared  to  K.p.  infection  ([13]  and  unpublished  observa¬ 
tion).  On  POD6,  after  antibiotic  killed  the  majority  of 
planktonic  bacteria  and  selected  for  biofilm,  P.a. -infected 
wounds  showed  a  network  of  inflammation/infection-re¬ 
sponse  molecules  predicted  to  be  regulated  by  type  I  in¬ 
terferons,  which  was  not  associated  with  enrichments 
indicative  of  a  stress  response  of  cutaneous  cells.  The 
network  of  inflammation/infection-response  molecules 
on  POD6  was  similar  to  a  network  unique  to  P.a. -in¬ 
fected  wounds  on  POD4  that  was  minor  part  of  the  ex¬ 
tensive  inflammation/infection-response  on  that  day.  On 
POD  12,  P.a. -infected  wounds  showed  indications  of 
greater  inflammation  and  impaired  healing,  which  may 
be  partly  related  to  a  larger  P.a.  biofilm  load  at  this  time 
as  compared  to  K.p.  biofilm,  as  suggested  by  our  previ¬ 
ous  studies  [13,49]. 

POD4 

Inflammation/infection-response 

Seeking  potential  mechanisms  of  the  very  similar  but 
much  more  extensive  inflammatory  host-response  of 
P.a.  vs.  K.p.  infection,  we  analyzed  predicted  upstream 
activating  molecules.  Those  from  bacteria  were  similar 


Table  3  Topmost  five  upstream  regulators  predicted  for  the  55  Upr  l/U-DEGs  unique  to  P.a.-infected  wounds  on 
POD6  in  the  categories  Cytokine,  G-protein-coupled  receptor.  Growth  factor.  Ligand-dependent  nuclear  receptor, 
and  Transmembrane  receptor 


Upstream 

Log 

Activation 

p  value  of 

Target  molecules 

regulator 

ratio 

z  score 

overlap 

in  dataset 

TLR3 

0.340 

3.111 

7.1 3E  19 

CXCL10,  DDX58,  DHX58,  GBP4,  GCH1,  IFIH1,  IFIT1,  IFIT3, 

IL15,  IL23A,  IRF7,  STAT1,  TLR3,  TNFSF10,  USP18,  ZNFX1 

IFN  Beta 

3.395 

1.34E  16 

CXCL10,  IFIH1,  IE  ITT  IFIT1B,  IL15,  IRF7,  STAT1,  TAPBP, 

TLR3,  TNFSF10,  USP18,  XAF1 

Interferon  alpha 

3.913 

1.67E  16 

CXCL10,  DDX58,  DHX58,  IFIH1,  IFIT1,  IFIT1B,  IFIT3,  IL15, 

IRF7,  PRF1,  STAT1,  TAPBP,  TLR3,  TNFSF10,  USP18,  WARS 

IFN  type  1 

2.087 

3.44E  15 

CXCL10,  DDX58,  DHX58,  IFIH1,  IFIT1,  IFIT1B,  STAT1,  TNFSF10,  UBA7 

Ifn 

2.793 

4.54E  12 

CXCL10,  DDX58,  DHX58,  IFIH1,  IFIT1,  IL15,  IRF7,  STAT1,  TLR3 
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for  both  infected  wound  types.  To  further  search  for 
specific  gene  expression  patterns  that  discriminate  the 
two  infected  wound  types,  we  analyzed  other  upstream 
regulators  in  categories  of  cytokines,  growth  factors  and 
signaling  receptors.  Twenty-five  upstream  regulators, 
predominated  by  type  I  interferons,  were  inferred  to  regu¬ 
late  56  I/U-DEGs  unique  to  P.a. -infected  wounds  that 
enriched  into  the  functions  Immune  Cell  Trafficking’ 
(36  molecules)  and  ‘Infection  (39  molecules;  ^-values 
1.24E-26  -  1.04E-7).  In  the  ‘Infection  category  twice  as 
many,  34,  enriched  into  the  subcategory  ‘Viral  Infection’ 
(p-value  2.56E-18)  as  enriched  into  ‘Bacterial  Infection 
(17  molecules,  p- value  8.37E-13).  This  network  could 
contribute  to  upregulated  IKBKE,  which  was  upregulated 
10-fold  in  P.a, -infected  wounds  but  only  2-fold  in  K.p. -in¬ 
fected  wounds.  IKBKE  is  a  non-canonical  IkappB  kinase 
with  an  important  role  in  regulating  antiviral  signaling 
pathways.  It  is  activated  downstream  of  the  recognition  of 
LPS  and  nucleic  acids  by  their  PPRs— TLR3  (dsRNA), 
TLR4  (LPS),  IFIH1  (long  dsRNA),  DDX58  (5’ppp  ssRNA, 
short  dsRNA),  and  ZBP1  (dsDNA)  [50]  —  and  activated 
IKBKE  phosphorylates  p65/RelA,  interferon  response  fac¬ 
tors  3  and  7  (IRF3  and  IRF7),  and  STAT1.  Consistent  with 
a  role  in  activating  IKBKE,  the  network  unique  to  P.a. -in¬ 
fected  wounds  on  POD4  includes  upregulated  TLR3, 
IFIH1,  DDX58,  IRF3  and  STAT1.  The  network  unique  to 
P. a. -infected  wounds  on  POD4  and  the  similar  network 
on  POD6  could  possibly  be  related  to  the  intracellular  rep¬ 
lication  of  P.a.  [46]. 

Some  gene  expression  induced  by  P. a. -infection  as 
compared  to  /Qz -infection  probably  involves  known 
virulence  mechanisms  of  P.a.  [46,51],  including  the  type 
III  secretion  effectors  ExoS  and  ExoT  that  are  expressed 
by  the  P.a.  strain  PAOl  that  we  used  for  our  studies 
[52].  A  previous  study  identified  genes  with  induced  ex¬ 
pression  dependent  on  ExoT  activity  [51]  that  are  some 
of  the  same  molecules  that  were  highly  expressed  in 
P.a.-  as  compared  to  TQz -infected  wounds  (e.g.,  ADM, 
DUSP5  and  NR4A1).  Additionally,  many  molecules  tar¬ 
geted  by  ExoS  activities  (Ras,  Ral,  Rabs,  and  ERM  pro¬ 
teins;  as  referenced  in  [53])  were  differentially  expressed 
in  P.a.-  vs.  K.p. -infected  wounds.  For  example,  the  ERM- 
like  molecule  NF2  was  one  of  the  36  genes  more  upreg¬ 
ulated  in  K.p.'  than  P.a. -infected  wounds. 

Stress-response 

In  parallel  to  the  greater  number  of  Upr-transcripts  and 
their  higher  upregulation  in  P. a. -infected  wounds  on 
POD4,  P.a. -infected  wounds  downregulated  more  tran¬ 
scripts,  and  most  of  the  common  Dnr-transcripts  were 
downregulated  to  lower  levels  as  compared  to  K.p.  -in¬ 
fected  wounds.  Dnr-transcripts  come  from  cutaneous 
cells,  since  cells  that  infiltrate  the  wound  in  response  to  in¬ 
fection  only  contribute  transcripts.  The  Dnr-transcripts  of 


infected  wounds  enriched  in  stress-response  pathways 
and  networks;  these  were  several-fold  more  extensive  in 
P.a. -infected  wounds. 

Given  that  a  constant  quantity  of  RNA  was  hybridized 
to  each  microarray,  more  leukocyte  infiltration  in  P.a .- 
infected  wounds  would  result  in  more  leukocyte  RNA 
and  proportionally  less  fibroblast  and  keratinocyte  RNA. 
But  this  does  not  appear  to  fully  explain  the  greater 
extent  of  Dnr-  transcripts  in  P.a.  wounds  on  POD4 
because  P.a. -infected  wounds  had  4-fold  more  Upr- 
transcripts  and  twice  as  many,  8-fold,  more  Dnr-transcripts 
than  K.p. -infected  wounds  (Table  1).  And  26%  of  overlap¬ 
ping  Dnr  I/U-DEGs  were  more  downregulated  in  K.p.-  than 
P. a.  -infected  wounds.  Also,  stress  responses  are  known  to 
involve  transcriptional  and  posttranscriptional  processes  in¬ 
cluding  the  specific  stabilization  and  destabilization  of 
mRNA  as  well  as  selective  translation  and  transcription, 
and  depend  on  the  particular  type  of  stress  [54,55]. 

The  downregulated  transcripts  assembled  into  net¬ 
works  with  highly  connected  nodes  such  as  CAND1  and 
cullins,  which  are  involved  in  protein  degradation,  and 
the  COP9  signalosome  subunits.  The  COP9  signalosome 
functions  in  the  regulation  DNA  repair,  checkpoint  con¬ 
trol  and  chromatin  remodeling  [41].  DNA  damage  was  a 
top  activated  function  enrichment  of  Dnr  I/U-DEGs 
unique  to  P.a. -infected  wounds,  which  is  a  likely  conse¬ 
quence  of  reactive  oxygen  species  from  the  strong  im¬ 
mune  response  to  P.a.  infection.  The  topmost  predicted 
upstream  regulator  was  MYC.  MYC  acts  as  a  transcrip¬ 
tional  amplifier,  inducing  expression  of  all  genes  to 
which  it  is  bound  [56].  The  predicted  inhibited  state 
of  MYC  in  P.a. -infected  wounds  is  consistent  with  gen¬ 
eral  suppression  of  transcription.  Also  upstream,  PPAR 
(peroxisome  proliferator-activated  receptor)  signaling 
was  predicted  to  be  in  an  inhibited  state.  PPAR  signaling 
regulates  many  proteins  in  epidermis,  which  is  very  ac¬ 
tive  in  lipid  metabolism,  and  where  all  PPAR  and  LXR 
(liver  X  receptor)  isoforms  are  expressed,  activating 
keratinocyte  differentiation  and  the  formation  of  the  per¬ 
meability  barrier  [57].  The  mTOR  pathway  was  also  pre¬ 
dicted  to  be  in  an  inhibited  state.  mTOR  is  a  conserved 
checkpoint  protein  that  senses  and  integrates  growth 
factor  signals,  and  levels  of  nutrients,  oxygen  and  energy. 
Major  downstream  targets  of  mTOR  signaling  are  com¬ 
ponents  of  the  translation  machinery  [58].  In  contrast, 
RARgamma  signaling  was  predicted  to  be  activated,  and 
RXRalpha-RARgamma  heterodimers  are  necessary  for 
growth  arrest  [59].  The  Dnr-transcripts  of  P.a. -infected 
wounds  enriched  4-fold  more  into  the  cell  death  func¬ 
tion  annotation  than  did  an  equal  number  of  Dnr- 
transcripts  of  K.p.  -infected  wounds. 

Previous  studies  have  suggested  that  eIF2  signaling  is 
activated  during  bacterial  infections  in  vitro  [60-63].  In 
epithelial  cells,  intracellular  bacteria  have  been  observed 
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to  up-regulate  several  proinflammatory  cytokine  genes 
(TNFA,  IL1B  and  IL6)  and  chemokine  genes  (IL8  and 
CCL20)  (as  referenced  in  [63]  and  therein).  These  genes 
were  among  the  most  highly  expressed  genes  in  our 
P.a.- infected  wounds.  Also,  human  vaginal  epithelial 
cells  (HVECs)  exposed  to  Staphylococcus  aureus  strongly 
upregulated  chemokines  and  cytokines  [64].  Among 
these,  on  our  microarray,  CCL20  was  significantly  upreg- 
ulated  only  in  P.a. -infected  wounds  (5-fold)  and  CXCL2 
was  upregulated  3-  and  14-fold,  respectively,  in  K.p.-  and 
P.a. -infected  wounds.  However,  our  data  cannot  discrim¬ 
inate  the  source  of  these  Upr-transcripts,  leukocytes  vs. 
cutaneous  cells. 

EIF2  signaling  was  among  the  top  three  canonical 
pathway  enrichments  for  Dnr  I/U-DEGs  of  both  infected 
wounds.  EIF2  is  phosphorylated  by  several  kinases  in  re¬ 
sponse  to  protein  disequilibrium  in  the  endoplasmic 
reticulum  (ER)  and  cytoplasm,  caused  by  diverse  cellular 
stresses  such  as  nutrient  deprivation,  exposure  to  the 
mTOR  inhibitor  rapamycin,  high  salt,  redox  perturb¬ 
ation,  and/or  DNA  damage  [65,66].  EIF2  phosphoryl¬ 
ation  shuts  off  global  translation  but  allows  synthesis  of 
specific  transcription  factors  and  their  target  genes  that 
support  recovery  from  the  stress  (e.g.,  ER  chaperones, 
oxidoreductases,  and  ERAD  components).  The  central 
transcription  factors  in  EIF2  signaling,  ATF4  and  ATF5, 
were  not  on  our  microarray,  but  the  related  transcrip¬ 
tion  factors,  ATF7IP,  ATF6,  ATF1,  BATF2,  ATF3,  ATF2, 
and  ATF7IP2  were  upregulated  in  P.a.-  relative  to  K.p .- 
infected  wounds. 

Beyond  a  threshold  of  stress,  apoptosis  can  result.  An 
apoptotic  switch  involving  the  Unfolded  Protein  Re¬ 
sponse  (UPR)  has  been  described  [67].  The  switch  to  the 
“terminal  UPR”  is  thought  to  involve  transcription  of 
the  pro-apoptotic  transcription  factor  CHOP  (DDIT3), 
which  was  upregulated  4.6-fold  in  7Q?. -infected  wounds 
and  13-fold  in  P. a. -infected  wounds  relative  to  unin¬ 
fected  wounds.  CHOP  has  been  suggested  to  induce 
death  by  promoting  protein  synthesis  through  increased 
transcription  of  GADD34  (PPP1R15A),  which  dephos- 
phorylates  EIF2  to  reverse  the  inhibition  of  translation, 
exacerbating  the  preexisting  ER  stress  [68,69].  Notably, 
in  our  study,  GADD34  (PPP1R15A)  was  30-fold  upregu¬ 
lated  in  P.a. -infected  wounds,  relative  to  uninfected 
wounds,  and  was  7.5-fold  more  highly  expressed  in 
P.a.-  than  K.p. -infected  wounds. 

POD6 

On  POD6,  the  wounds  containing  biofilm  expressed 
magnitudes  fewer  inflammation  transcripts  than  before 
antibiotic  treatment.  Although,  10-fold  fewer  P.a.  than 
K.p.  CFU  were  recovered  from  wounds  on  POD6  [12], 
only  P.a. -infected  wounds  contained  Upr  I/U-DEGs 
that  enriched  in  inflammation/infection  functions.  This 


network  was  virtually  the  only  network  found  on  POD6. 
It  was  similar  to  a  network  on  POD4  consisting  of  Upr 
I/U-DEGs  unique  to  P.a. -infected  wounds,  which  was  a 
minor  part  of  the  extensive  inflammation/infection-re¬ 
sponse  on  POD4.  In  common  between  the  two  networks 
were  three  identical  molecules  and  three  highly  similar 
molecules,  and  both  networks  contained  molecules  in¬ 
volved  in  nucleic  acid  sensing  with  inferred  regulation 
by  type  I  interferons.  However,  type  I  interferons  were 
not  upregulated  on  our  microarrays.  One  possible  ex¬ 
planation  for  these  networks  is  that  they  are  a  conse¬ 
quence  of  growth  of  P.a.  within  epithelial  cells  [46]  and 
could  be  regulated  via  an  interferon-independent  mech¬ 
anism  such  as  found  in  fibroblasts  [70,71].  Another 
possibility  is  that  these  transcripts  could  come  from  a 
specific  subpopulation  of  infiltrating  cells  such  as  den¬ 
dritic  cells  [72].  Consistent  with  the  greatly  reduced  in¬ 
flammation  on  POD6,  the  Dnr  I/U-DEGs  of  cutaneous 
cells  did  not  significantly  enrich  into  annotated  func¬ 
tions  suggesting  a  stress  response. 

Following  two  days  of  antibiotic  treatment,  POD6  wounds 
differentially  expressed  fewer  genes  than  on  POD4  or 
POD  12.  However,  a  large  fraction  of  the  few  I/U-DEGs 
present  in  P.a. -infected  wounds  were  related  to  each  other 
through  regulatory  connections,  including  inferred  upstream 
regulation  by  type  I  interferons.  Such  a  network  of  tran¬ 
scripts  could  possibly  be  an  identifying  marker  of  P.a.  bio¬ 
film  colonization  that  has  potential  to  develop  into  chronic 
infection  that  impairs  healing. 

POD12 

Relative  to  uninfected  wounds,  both  types  of  infected 
wounds  on  POD  12  differentially  expressed  genes  that 
enriched  into  the  functions  ‘Vasculogenesis’  and  Angio¬ 
genesis’,  with  predicted  activated  states,  and  ‘Connective 
Tissue  Disorders’.  These  similar  functions  between  the 
two  types  of  infected  wounds  are  likely  related  to  both 
being  open  as  compared  to  the  uninfected  wounds  that 
were  closed  on  POD  12. 

Transcripts  that  were  more  highly  expressed  in  P.a.~ 
as  compared  to  K.p. -infected  wounds  enriched  in  func¬ 
tions  and  networks  related  to  the  inflammatory  response 
and  connective  tissue  disorders,  suggesting  that  P.a. -in¬ 
fected  wound  were  more  inflamed  and  disordered. 
Additionally,  FGF1  (fibroblast  growth  factor  1)  was  the 
topmost  upstream  regulator  of  these  molecules  more 
highly  expressed  in  P.a. -infected  wounds.  FGF  family 
members  possess  broad  mitogenic  and  cell  survival  ac¬ 
tivities,  and  FGF1  is  a  known  modifier  of  endothelial  cell 
migration  and  proliferation,  as  well  as  an  angiogenic 
factor.  However,  as  an  upstream  regulator  of  these 
DEGs,  FGF1  was  suggested  to  be  in  an  inhibited  state. 
Additionally,  PLI<5  (polo-like  kinase  5)  was  upregulated 
10-fold  in  P.a.-  relative  to  K.p. -infected  wounds.  PLI<5 
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has  been  shown  to  induced  G1  arrest  when  overex¬ 
pressed,  to  be  downregulated  in  proliferating  cells,  and 
to  accumulate  in  quiescent  serum-starved  fibroblast 
[73].  These  results  could  indicate  an  impaired  state  of 
healing  of  P.a. -infected  wounds.  Conversely,  upstream 
regulators  of  the  more  highly  expressed  transcripts  in 
K.p.-  vs.  P.a, -infected  wounds  included  a  glucocorticoid 
and  PDGF.  These  predicted  upstream  regulators  may  in¬ 
dicate  lower  inflammation  (glucocorticoid)  and  active 
healing  (PDGF)  [74].  Overall  these  expression  data  and 
their  regulator  and  function  enrichments  suggest  that 
P.a. -infected  wounds  had  greater  inflammation  and  K.p.- 
infected  wounds  had  more  developed  and  active  healing. 

Conclusions 

These  data  contribute  to  defining  the  immune  response  of 
wound  and  proximal  skin  to  bacterial  infections  of  different 
virulence.  In  the  clinical  setting,  diagnosis  of  infected 
wounds  can  be  complicated  by  an  absence  of  signs  of  in¬ 
flammation  (e.g.,  as  on  POD6)  and  late  diagnosis  can  result 
in  impaired  healing  [75,76].  Current  diagnosis  of  infection  is 
based  on  recoverable  CFU  which  captures  only  culturable 
bacteria.  PCR  identifies  more  bacteria  in  wound  samples, 
but  mixed-species  infections  could  have  unpredictable  viru¬ 
lence.  The  host-response  counterbalances  infections,  and 
molecular  markers  of  the  host-response  to  infection  may  ul¬ 
timately  be  useful  for  diagnosis  and  treatment.  The  current 
study  suggests  that  patterns  and  extents  of  expression  of 
Dnr-transcripts,  expressed  in  cells  resident  to  the  wound 
and  proximal  skin,  as  well  as  specific  Upr-transcripts,  may 
contribute  to  the  prediction  or  identification  of  poorly 
healing  wounds,  and  may  be  useful  for  indicating  wound- 
specific  treatments  to  improve  wound  outcome. 
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